Dynamics of Two Qubits: Decoherence and an Entanglement Optimization Protocol 
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The evolution of two qubits coupled by a general nonlocal interaction is studied in two distinct 
regimes. In the first regime the purity of the individual qubits is interchanged through the entangle- 
ment shared by the two. We illustrate how this can be a mechanism for decoherence. In the second 
regime, the interaction entangles two initially pure qubits. The dynamical maps for the reduced 
unitary evolution of both initially simply separable and not-simply-separable states are found. We 
outline a protocol for optimizing the entanglement generation subject to constraints. 



I. INTRODUCTION 

Universal, two qubit quantum gates are the basic units 
from which quantum information processing devices may 
be constructed [l], Physical implementations of two 
qubit gates depend on understanding and controlling 
the possible interactions between the individual qubits. 
Zhang et. al. [3j have obtained the most general form 
for non local two-qubit interactions. A complete picture 
of the dynamical possibilities allowed by this interaction 
is yet to be pieced together, although some numerical 
efforts have started to scratch the surface [1]. Entan- 
glement is a resource for a quantum information proces- 
sor [j| Q, while entanglement of elements of the infor- 
mation processor with its environment leads to undesir- 
able loss of purity and decoherence. Naturally, there has 
been a lot of interest 1 in understanding M , quantifying 
@, ©, El El El, and controlling [H [ll 111 entangle- 
ment and also preventing decoherence. We show that 
both the loss of purity of qubits and the generation of 
entanglement between two qubits in an optimal way can 
be understood in terms of the dynamics generated by the 
same general two qubit interaction under different con- 
ditions. 

Our main result consists of an explicit expression for 
the dynamics of two qubits coupled by the general Hamil- 
tonian. We calculate the dynamical map for both the 
initially simply separable and not simply separable two- 
qubit state. We focus on how the dynamics affect the 
purity of the individual qubits. We describe purity swap- 
ping, where one qubit is dynamically purified at the ex- 
pense of another. This may be contrasted with entan- 
glement swapping We also study the generation 
of entanglement by the dynamics. In [lj} a method for 
maximizing the rate of entanglement generation has been 
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1 "As the strong man exults in his physical ability, delighting in 
such exercises as call his muscles into action, so glories the analyst 
in that moral activity which disentangles." [sic] Poe, E.A. "The 
Murders in the Rue Morgue" Complete Stories and Poems of 
Edgar Allan Poe. Ed. Doubleday & Company, Inc. Garden 
City, NY (1966) 



proposed given some entangling Hamiltonian, and the 
capability to do instantaneous local operations on each 
qubit. We propose a protocol that optimizes any given 
entangling process subject to reasonable constraints that 
may be found in experiments. 

The most general dynamics of the quantum state rep- 
resented by a density matrix p can be described in terms 
of a dynamical map [H, [l9[ : 



B(p) 



(1) 



The dynamical map represents the effect of the coupled 
unitary evolution of the system and its environment [201 ] . 
In other words: 



P 



p A (t) = Tr B [Up J 



P 



(2) 



where p A is the state of the system and p B is that of the 
environment. In the case where the system and the en- 
vironment are finite dimensional and start in a separable 
state, p A <£> p B , 23 is a matrix with positive eigenvalues 
and it is said to be completely positive. 



II. EVOLUTION FOR TWO INITIALLY 
SIMPLY-SEPARABLE QUBITS 

We concentrate on the simplest example, where both 
p A and p B are single qubits: 

,A , „.„A ,B ^ b . a B 



P A = 



1- 



A «± 
P B = ~ 



(3) 
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with <j A and of being the Pauli spin matrices for each 
of them. The Bloch vectors a = (ai, 02,03) and b = 
(61,62,^3) provide a convenient way of parameterizing 
single qubit states that we will use in this paper. To- 
gether, p and p B form the initially separable 4x4 state, 



E AB (0) = p A ®p B = ~ [t J 



t B + hl A ® a? ) 



Y aibjof (g> a B \ 

i,3 



(4) 



where subscripts take values from {1,2,3}. The most 
general Hamiltonian for two qubits is: 



H 



a 



2 



In the interaction picture, it becomes, 



hi 



which has nine parameters. Using local unitary trans- 
formations with three parameters each, the number of 
parameters can be brought down to three [H, S3, SI SI: 



H(t) -» H(t) = 5>i(*)<tf ® at 



(7) 



The time evolution of the overall state E AB is given 
by, E AB {t) = UE AB (0)U\ where U = Te" 4 J ' H ^ dt . For 
simplicity, we will assume that there is no free evolution 
for individual qubits, making 7,(i) — » 7i, and, 

3 

U = II [ cos (7^)l A ® 1 S - i sm{jjt)af <g> of] . (8) 

3=1 



To calculate E AB (t), we use the property that oy (g> 
i 4 , (gi o B , (X) cr^ all commute with each other. For 



(6) each of the terms of E AB (0) we obtain: 



Uoi ® 

Ul ® O.U ] 

Uoi (g> a j [/ f 



10 1, 

<J; (g) O; 



= 1 



cr, : e 



,2«(7(. 



5<Ti +7j<Tj 



(9) 



Using Eq. (|9|) , the evolution generated by Eq. 
two qubit density matrix is: 



of the 



E AB (t) = - [l A ® 1 S + (Qftaf (8 I s + S,£*1 A (8 af + C k S 3 o A ® of - CjS k of <g> erf) 

i=i 

+6 s: (QCfel- 4 O erf + SjSfcO^* ® I s + (8 of - CjSfcO^ <g> of) + a^of <g> of 

+ ai bj (CiCjof <g> of + SiSjof ® of + C^o-fc <g> I s - C 3 Sil A ® of ) 

+a j b i (C i C j a* ® erf + S^-of ® of + C^l^ 4 ® of - C^of ® I s )] , (10) 



where C; = cos (2t7i) , 5* = sin (2*7, ) and k are 
cyclic. In other words, the coefficients of the Pauli ma- 
trices in Eq. JU) transform as follows: 



h 



cik- 



Si Sj 



Qjbi i f< q 
a,h, ' ' ; S; 



(11) 



where {i,j, k} are cyclic, and 



a-ibj 



a i bjC l Cj + ajbiSiSj 
+Cijk(bkCjSi — akCiSj), 



(12) 



where it is not required for {i, j, k} to be cyclic or distinct. 



III. REDUCED DYNAMICS OF TWO 
INITIALLY SIMPLY SEPARABLE QUBITS 

To find the reduced dynamics of the system, p A , we 
just need to carry out the partial trace from Eq. ||5J). Us- 
ing Eqs. |T8]). (fPTj) . ([12]) and the fact that of are trace- 
less, this is quite straightforward: 



p A (t) 



l + a 3 (t) ai(i) - ia 2 (t) 
ai(t) + ia%{t) 1 — a 3 (t) 



(13) 



where 



a,i (t) = atCjCk + biSjS k 

+ djbkCjSk — akbjCkSj. 



(14) 



The purity of the reduced density matrix changes with 
time: 



P A {t) 



Tr 

1 r 

2 



p A W 



1 + a x (tf + a 2 (t) 2 + a 3 (t)' 



(15) 



Following a similar procedure we can find the evolution 
of p B (t) and its purity P B (t). The time evolution of p A 
can be described using the dynamical map from Eq. (fT]). 
The mapping matrix can be explicitly calculated to be 
El: 



3 



$g(o,t) 



/ 1 + b 3 S 1 S 2 + dC 2 b 2 C 1 S 2 -ib 1 C 2 S 1 (hSs - b 2 C 3 ) S 2 (Ci + C 2 ) \ 

-i(b 2 S 3 + b 1 C 3 )S 1 x(C 3 +ib 3 S 3 ) 

b 2 CxS 3 + ib l C 2 S 1 1 + 6 3 5i5 2 -CiC 2 (C 2 - C x ) (biS 3 + b 2 C 3 ) S 2 

x(C 3 ~ib 3 S 3 ) -i{b 2 S 3 -b l C 3 )S l 

{b 1 S 3 -b 2 C 3 )S 2 (C2-C1) 

+i(b 2 S 3 + b 1 C 3 )S 1 x(C 3 + ib 3 S 3 ) 1 - 6 3 SiS 2 - C X C 2 -b 2 C x S 2 + ib x C 2 S x 
(C1 + C2) (6i5 3 + 6 2 C 3 )5 2 

V x (C 3 - i&a-Sa) +i (6 2 S 3 - &1C3) S ± -b 2 dS 2 - s&iC 2 Si 1 - 6 3 5iS 2 + CiC 2 / 



(16) 



Note how the dynamical map carries the influence of ply separable, i.e.: 
p B through the parameters b\, b 2 , b 3 . Thus, given any 
known initial state p A and sufficiently detailed evolution 
p A (t), its interaction with another unknown p B can be 
reconstructed, and even used to determine the parame- 
ters for the unknown state. 

The most general map on a qubit can be implemented ^ 1 . 

as the contraction of the unitary evolution of the given — 4 2^/ L + e 4 ocr 4 ® + 

qubit coupled to at most two other qubits. We have re- 1,0 

stricted to the case where there is only one other qubit, e oil ® o» + ® °i] 1 (17) 

thereby excluding certain maps. If p B is allowed to be 
a mixed state, the family of dynamical maps we are ex- 
cluding by choosing a one qubit environment is very small 



IV. DYNAMICAL MAP OF THE GENERAL 
TWO QUBIT-STATE 

Eqs. (|11I12[) may be extended in a straightforward 
fashion to the case where the initial state, E, is not sim- 



where it is not necessary that = e^o x eoj. This 
would include initially entangled states. The evolution 
for the state E under Eq. (jSj) can be computed using the 
result from Eq. ([9]): 



E{t) = - [l A ® ! S + e»o {CjCkO-f ® I s + S^l^ 1 erf + C fc Vfc ® of - C^af ® of) 

+e ; (QCfel^ 4 <g) erf + SjSkCrf ® 1 S + C k S-jaf ® erf - CjS k a A <g> erf) + e«o^ ® erf 
+e y - (CiCjcrf <g> af + S^crf <g> erf + C l( S>fc ® I s - Cj^l- 4 ® erf) 

•( ,,((",(>/ ® af + Si^oi 1 ® af + CiSjl^ (g) fff - CjSitT^ ® 1 B )] . (18) 



Again, to find the reduced dynamics of the system we 
just need to carry out the partial trace, and the evolution 
of each of the components become: 

eio (t) = eioCjCk + eoiSjSk 

+ ejkCjSk - ekjCkSj. (19) 



As before, we would like to construct the dynamical 
map for this evolution. This time, the map has to carry 
the parameters of the trace out qubit as well as the cross- 
terms. Some of the terms that could be factored before 
into parameters in the reduced state and parameters in 
the map cannot be. This leads to the following map: 
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1 + CiC 2 

+60351 ^2 

■\-GviG\Si — e2iC , 2iS'i 


£015253 + ie 25 3 S , i 
+e23C253 — e 3 2C 3 S2 
+ie 3 iC 3 Si - ie 13 dS 3 

{C X +C 2 )C 3 





l-CiC 2 
+eo3Si<S , 2 

■\-e\2C\S2 — &2\C2S\ 

(C2 — C\) C 3 

6015253 + ieo25 3 Si 
+e23C25 3 — e32C 3 S , 2 
+ie a iC 3 Si - iei 3 CiS 3 



eoi5 2 5 3 - ie 2 5 3 S'i 
+623^53 — ez2C 3 S 2 
-izziCzSi + ie 13 CiS 3 



(d + C 2 ) C 3 



eoi5 2 5 3 - ze 25 3 S'i 
(C2 — Ci)C3 +623^53 — e 3 2C 3 S2 

ie 3 iC 3 Si + iei 3 C\S 3 



1 - CiC 2 
— eo35i52 
-ei 2 C\S2 + e2iC25i 





1 + CiC 2 
— 603516*2 
-e^CiS^ + e 2 iC 2 5i / 



The first thing to note about this map is that there 
are some elements that are 0. This is due to the fact 
that not all terms can be decomposed as a product and 
split between the state and the action of the map. In 
other words, the map 25 carries all the information in 
the cross terms of the initial two qubit state. To guaran- 
tee that this map correspond to something physical, we 
need to observe that due to the initial correlations of the 
bipartite state, only certain values are permitted for the 
reduced intimal state. The map has some information of 
its allowed domain, and as long as it acts on it the evolu- 
tion can still be given physical interpretation [241 ] and its 
eigenvalues in general can be negative. There are some 
experimental examples of these non-completely positive 
maps [25] and also their role in decoherence control tech- 
niques [261 ]. 



Purity 




4 2 4 



FIG. 1: Purity as a function of time for a\ = 72 = 73 = 
1, 02 = 03 = b\ = 62 = &3 = 7i = 0. The solid line represents 
P A (t) while the dashed line represents P B (t). At t = 7r/4 
purity has been totally swapped. 



V. DECOHERENCE MODELLING 

Now, we go back to the case of initially simply- 
separable states. We want to study certain kinds of evo- 
lution to show the use of these maps in some decoherence 
modelling schemes. 

Consider the case where p A (0) is pure while p B (0) is 
fully mixed. Assume some interaction where the only 
non-zero parameters are a\ = 71 = 72 = 73 = 1. This 
interaction is chosen to be of the "swap-gate" form, but 
many interactions can have a similar effect. Figure [T] 
shows how the purity of p A (t) and the purity of p B {t) 
change with time. As time progresses, they start to get 
entangled and exchange purity through entanglement. 
P A (t) and P B (t) are equal at t — n/8, where some purity 
has been lost to entanglement. At t = ir/4 they become 
separable again (2?| purifying p B at the expense of p A , 
a dynamical process we call purity swapping. If only a 
certain component of the qubit is measured and no cor- 
relations studied, this might look like Rabi oscillations, 
although the Bloch vector is not rotating, but oscillating 
in length. 

We remark that for a weak coupling of this type, at 



very short times, P A (5t) can only decrease (except in 
trivial cases). This can be considered to be a mechanism 
of decoherence. We can model a reservoir as a stream 
of {p B }, where each of them interact independently for 
a short average time [28j . swapping some purity from 
p A to each pf , but stopping the coupling before there is 
enough time to return the stolen purity. We can think 
of this as a quantum version of the Boltzmann gas [2!| . 
This corresponds to acting with the dynamical map from 
Eq. ((THJ) in sequence: 

p A S (t„-i,t„) f B (i„- 2 ,t„-i) . .. Q3(W0 (^) . (21) 

By controlling the strength, duration and number of 
these reservoir interactions it is possible to model deco- 
herence processes, as well as desired, using only a finite 
number of degrees of freedom for the reservoir. 

There are many cases where the periodic couplings are 
a good approximation and its properties exploited to pre- 
vent decoherence. Bang-bang control, for example, is a 
technique where the periodicity (on average) of a deco- 
herence coupling can be synchronized with a local control 
pulse effectively creating a decoherence-free space for the 
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state of interest [3(1 [3l| . 

VI. OPTIMAL ENTANGLEMENT 
GENERATION 

Another interesting regime to study is related to the 
creation of maximally entangled Bell states. Assume that 
we have initially two pure states, {cii = 1, a 2 = a 3 — 
0}, {62 = 1, bi = 63 = 0}. At a specific time tbeii, each 
of the qubits' purity goes to a minimum. For 73 = 1, the 
minimum is at tbeii = tt/4, and we get that E AB (tbell) is: 

- [t A ®t B + a A ®a%-<j£®a% + a£® crf] , (22) 

that, given some freedom to choose the basis for p B , 
would be equivalent to the Bell state (|00) + |11)). A 
similar, experimentally feasible time-reversed procedure 
would be responsible for extracting purity out of entan- 
glement (32|. 

How much can an entanglement creation experiment 
be improved? From Eq. (Tj"5|) an entanglement optimiza- 
tion protocol can be determined, given certain reasonable 
constraints. In particular, we would like to show how 
given a certain two-qubit state and certain experimen- 
tal limitations in the Hamiltonian that can be applied, 
the evolution can be optimized to maximize the entan- 
glement between them. 

Lets illustrate this with a simple example for the case of 
the creation of entanglement from initially pure qubits. 
Limitations on the allowed interaction and their dura- 
tion would come from particular experimental require- 
ments. This example, although very simple, contains all 
the elements to illustrate how a procedure like this can 
be implemented. The more general protocol is discussed 
afterwards. 

Imagine some experimental setup that prepares two 
qubits with the Bloch vectors a and b representing their 
states oriented at some angle with respect to each other. 
The aim is to entangle these two. Assume that only 
two kinds of couplings between them are allowed and we 
can only act first with one kind followed by the other, 
due to some experimental constraints. The only thing 
that can be controlled is the duration for which each 
coupling is used. The total time available is restricted 
by decoherence. This type of interaction can be similar 
to the "strongly modulating" pulses in nuclear magnetic 
resonance implementation of quantum algorithms, but is 
not limited to it. How long should we act with each of 
them to maximize the entanglement rate? 

Let the two qubits p A and p B , initially be in the states 
given by a\ — 1, 02 = 03 = and b\ = b% = 1/ V2, 63 = 0. 
The allowed interactions are 72 = 1, 71 = 73 = for 
some time [0, t'] and followed by 73 = 1, 71 = 72 = at t' 
for an interval (t! ', n] . Using Eq. I|18p we can calculate the 
state E AB (t) of the system at t' . Using E AB (t') as the 
new initial condition in Eq. I|17p and the new coupling, 
73 = 1 we can compute the state of the two qubits during 



(£', 7r] . Eq. (p~3|) gives us the reduced density matrix of 
one of the qubits as a function of time from which we can 
compute its purity at every time. 

What is the time t' that gives us the maximum en- 
tanglement? This protocol is not dependent on a par- 
ticular measure of entanglement, but we will choose for 
simplicity the entropy of entanglement ||, which is a 
good measure as long as E AB (t) remains pure. Since it 
is monotonically related to the linear entropy [13, l33l |. 
the entanglement measure £ can be also chosen: 

£ ~ 1 - Tr [(/) 2 ] = 1 - Tr [{p B f} = 1 - P. (23) 

Using Eqs. (fl4|) and (fl~5|) . we find that for our choices of 
the parameters: 

^^ 3003(8^ + 29 (24) 
32 

which is maximized at t' = 7r/8, 37r/8, 57r/8 or 77r/8. 
We admit that this example is rather simple, and was 
chosen because it can be easily solved algebraically. More 
complicated and realistic couplings and limitations can 
be solved numerically with similar results. 

This method can be applied to any interaction of the 
form Eq. ([7]) , allowing us to optimize a very general class 
of entanglement creation procedures. It is significantly 
different from the procedure proposed in [l7| , in that we 
do not need to assume full control over the local trans- 
formations on each qubit at all times. 

The general scheme of this entanglement optimization 
protocol can be summarized as follows: 

1. Choose the initial conditions for the system, and 
express them in the form of cither Eq. (0| or 

Eq. jnj. 

2. Identify the coupling parameters and constraints 
for Eq. ([7]). Compute the evolution due to the uni- 
tary operator from Eq. ([5]) by using the transfor- 
mations in Eqs. (|llll2U8p . 

3. Choose the measure of entanglement of your pref- 
erence @, i m m d . 

4. Considering the experimental constraints, optimize 
with respect to the desired parameters. 

VII. CONCLUSIONS 

In summary, we have explicitly calculated the evo- 
lution generated by a general interaction between two 
qubits, opening the doors for studying all sorts of entan- 
gling interactions and experimentally realizable universal 
quantum gates. We compute the Dynamical Map of the 
most general evolution for two qubits, and showed the 
different nature that the Map can have when the state is 
initially simply-separable of not. Through examples, we 
illustrated how purity and entanglement are interchange- 
able quantities. We studied dynamical purity swapping, 
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and its connection to decoherence phenomena. Our sim- 
ple two qubit interaction is flexible enough to study the 
fundamentals of decoherence that are usually studied as 
an interaction with an infinite degrees of freedom reser- 
voir. We described a procedure to generate Bell States 
from pure states, and its converse, purification. Finally, 
we proposed a practical protocol to optimize a given en- 



tangling procedure under realistic constrains. 
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